Cluster aggregation model for discontinuous percolation transition 
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The evolution of the Erdos-Renyi (ER) network by adding edges is a basis model for irreversible 
kinetic aggregation phenomena. Such ER processes can be described by a rate equation for the 
evolution of the cluster-size distribution with the connection kernel Kij ~ ij, where ij is the product 
of the sizes of two merging clusters. Here we study that when the giant cluster is discouraged to 
develop by a sub-linear kernel Kij ~ (ij)" with < uj < 1/2, the percolation transition (PT) is 
discontinuous. Such discontinuous PT can occur even when the ER dynamics evolves from proper 
initial conditions. The obtained evolutionary properties of the simple model sheds light on the origin 
of the discontinuous PT in other non-equilibrium kinetic systems. 

PACS numbers: 64.60.ah,02.50.Ey,89.75.Hc 



Irreversible cluster aggregations are widespread phe- 
nomena occurring in a diverse range of fields, including 
dust and colloid formation, aerosol growth, droplet nuclc- 
ation and growth, gelation transition, etc [l| . The Smolu- 
chowski coagulation equation can successfully de- 

scribe such cluster aggregation processes. In linear poly- 
merization, molecules with two reactive ends can react 
to form long chains. In this case, the reaction kernel is 
given as Kij = 1 , where i and j are the masses of the two 
reactants. For the aggregation of branched polymers, the 
reaction kernel has the form Kij = (ai + b)(aj + b), where 
a and b are constants. When clusters have a compact 
shape, the reaction kernel has the form Kij ~ (ij) 1 ~ 1 ' d , 
where d is the spatial dimension. Intensive studies have 
been carried out using the Smoluchowski coagulation 
equation with such different kernel types [l|, UQ, and 
it is known that sol-gel transitions can occur at cither 
finite or infinite transition points. They are continuous 
transitions. 

During the past decade, the evolution of complex net- 
works has been of much interest to the science commu- 
nities in multidisciplinary fields. To study percolation 
transition (PT) during network evolution, the branch- 
ing pro cess approach jjj [l(| and the Potts model formal- 
ism have been used. Such complex network evolution 
can also be viewed as a cluster aggregation phenomenon, 
and can be studied by the rate-equation approach Q. 
For example, in the evolution of the classical random 
network, called the Erdos-Renyi (ER) model, an edge 
is added at each time step, thereby either connecting 
two separate clusters (inter-cluster edge) or increasing 
the edge number in one cluster without changing clus- 
ter numbers (intra-cluster edge). Fig. [T] shows that the 
frequency of inter-cluster connections is dominant until 
the percolation threshold. Thus, the cluster aggregation 
picture of the ER network evolution comes in naturally. 
In this paper, we extend the cluster aggregation dynam- 
ics in networks to more general cases. Specifically, the 
model is as follows: In a system composed of N vertices, 
we perform the following tasks at each time step. 

• Two clusters of sizes i and j are chosen with proba- 
bilities qi and qj , respectively. The two clusters can 



be the same. Probability qi is given as ki/ k s n s , 
where fcj and Ui are the weight and density of an 
i-sized cluster, respectively. 

• Two vertices are selected randomly one each from 
the selected clusters. If they are not yet connected, 
then they are connected by an edge. If they are 
already connected, we choose another pair of ver- 
tices in the same manner until a link can be added. 
Self-loop cases are excluded. 

We repeat these simple steps until a given time t = L/N, 
where L, the number of edges added to the system, is 
tuned. This model is called the cluster aggregation net- 
work model hereafter. In this model, the two selected 
clusters can be the same, and thus, the evolution can 
proceed even after one giant cluster remains. The ER 
network corresponds to the case ki = i. Here, we show 
that when the weight is sub-linear, as ki — i u with 
< oj < 1/2, a discontinuous PT occurs at a finite 
transition point. Moreover, under certain initial con- 
ditions, the ER dynamics also exhibits a discontinuous 
PT. This observation is remarkable, because a discontin- 
uous PT has rarely been discovered in irreversible kinetic 
systems, except for recent observations in the ER [l3| 
and other networks fl4l - fl9} under the so-called Achliop- 
tas process [l3[. On the other hand, it is noteworthy 
that the cluster aggregation network model evolves by 
single-edge dynamics, as compared with the ER network 
under the Achlioptas process, which involves a pair of 
edges at each time step. Thus, this cluster aggregation 
network model allows us to study the underlying mecha- 
nism of the discontinuous PT analytically for some cases, 
which is shown later. 

The cluster aggregation processes in the model are de- 
scribed via a rate equation for the cluster density, which 
takes the following form in the thermodynamic limit: 



dn s (t) 
dt 
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where c(t) = ^2 s k s n s (t). The connection kernel K^ = 
kikj/c 2 . The first term on the right hand side repre- 
sents the aggregation of two clusters of sizes i and j with 
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FIG. 1: (Color online) The fraction of each type of attached 
edges, inter-cluster (•) or intra-cluster (v) edges for the ER 
model. Arrow indicates percolation threshold at p c = 1/2. 



i + j = s and the second term represents a cluster of 
size s merging with another cluster of any size. The 
rate equation differs from the Smoluchowski coagulation 
equation in two aspects. First, the connection kernel is 
time-dependent through c(t) when u) ^ 1. Second, the 
second term on the right hand side of Eq. ([T]) includes the 
process of merging with an infinite-size cluster. Hence, 
Eq. ([1) with ui = 1 and c = 1 describes the ER pro- 
cess, while the conventional Smoluchowski coagulation 
equation with ui = 1 does not, because only sol-sol reac- 
tions are taken into account. However, the case including 
the infinite-size cluster in the Smoluchowski coagulation 
equation was also considered in Ref @ , which was called 
the F-model. Owing to the presence of c(t) , a PT occurs 
at a finite transition point even when a PT does not occur 
in the Smoluchowski coagulation equation, for example, 
when ui = 0. Here, we study the cases fcj = 1 (uj = 0), 
fcj = i u with < ui < 1, and fcj = i (u> = 1), separately. 

The case ui = 0: In this case, c{t) = ^ s n s becomes 
the total density of the clusters, which decreases linearly 
with time. The generating function of n s (t) is defined 
as f(z 7 t) — J2 s ns (t) zS > wner e z is the fugacity in the 
range < z < 1. Then, one can obtain the differential 
equation for f(z,t) from Eq. ([T]) and solve in a closed 
form as }{z,t) = (1 - t) 2 z/(l - zt) for t < 1 and for 
t > 1 in the thermodynamic limit. Expanding f(z,t) as 
a series in z, we obtain 



n s (t) = (l-t) 2 t 



2.S-1 



(2) 



for t < 1. This formula shows that the cluster size dis- 
tribution decays exponentially as s becomes large. Par- 
ticularly, when 5 = 1 — t is small, n s (8) « S 2 e~ s ^ s with 
s* w 1/(5. The characteristic size s* diverges as 5 — > 0. 
As shown in the inset of Fig. [2ja) , n s (t) is almost flat at 
S = 1CP 4 for = 10 6 , indicating that large-size clusters 
are relatively abundant. The merging of these clusters 
causes a sudden jump in the giant cluster size, leading to 
a first-order transition. 

We find the giant cluster size G(i) by using the relation, 
G{t) = l — f'(l,t) = 1 — J2' s sn s (t) , where the summation 
excludes an infinite-size cluster. We find that 



G(t) 



Off < t < 1, 
Iff t > 1, 



(3) 
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FIG. 2: (Color online) The cluster-size distributions of the 
cluster aggregation network models with (a) k% = 1, (b) hi = 



(c) ki 



, and (d) ki — 1 (ER model) with mono- 



disperse initial condition, (a) is drawn in semi-logarithmic 
scale, while the others are in double-logarithmic scales. Data 
points in (a) are at t/t c = 0.3(»), 0.5(V) and 0.7(o), and 
in (b), (c), and (d), they are obtained at t/t c = 0.50 (•), 
0.95 (V), 0.998 (o)), and 1.003 (o)). Data points in the inset 
of (a) are at <5 = 1 — t = 10 -4 . In (b), there exists a hump for 
data points (o). The system size iV is taken as 10 6 for (a) and 
10 for (b), (c), and (d). Solid lines in (b) and (c) represent 
analytic formulae Eq. (Q 



in the thermodynamic limit (Fig. [3Ja)). Thus, the PT 
is first-order at t c = 1. This result differs from what we 
obtain from the Smoluchowski coagulation equation, in 
which the transition point t c = oo. 

The case < ui < 1: For this case, while exact solu- 
tion for n s (t) is not obtained, n s (t c ) is done under cer- 
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FIG. 3: (Color online) The density of the giant cluster size G 
versus time t showing (a) discontinuous transitions in ki = 1 
(o), ki — i ' 2 (A), ki — i ' 4 (•), and (b) continuous transitions 
in the cluster aggregation network model ki = i 0,6 (V), ki = 
i°- s (A), and the ER network (o). N = 10 5 in both (a) and 
(b). 



tain assumptions. To proceed, we define the generat- 
ing function g u (fi,t) = J2 S sU ' n s{t)e IJ,s /c(t) ([i < 0), and 
presume that n s (tc) ~ s~ T . Next, we use the assump- 
tion made in Ref. |5j, [7| for the Smoluchowski coagulation 
equation that n s (t) = n s (t c )/(l + b(t — t c )) near t = i+, 
where b is an s-independent constant. Then, comparing 
the most singular terms in the series of the generating 
functions /(e^ic) and g^(/j,,t c ) in /i, we find that 



1 + 2uj 
3/2 + u. 



< w < 1/2, 
1/2 < u < 1. 



(4) 



This result is confirmed numerically in Fig [5] When t < 
t c , n s (t) follows a power-law function with an exponential 
cutoff for 1/2 < uj < 1, but it exhibits a hump in a large- 
size region for < u) < 1/2 (Fig. [^b)). 

We examine G(t) as a function of time for various to 
cases. G(t) exhibits a transition at finite t c , which is 
continuous for 1/2 < uj < 1, discontinuous for < u> < 
1/2 (Fig. [3]), and marginal for uj = 1/2. The first-order 
transition is tested in Fig. [4] using the scaling approach 
introduced in Ref. [llj]. We define A = <i — to, where to 
and t\ are chosen as the times at which the value of G(t) 
reaches 1/VN and 0.8 for the first time, respectively. We 
find numerically that for < uj < 0.5, A decays as TV — s- 
oo, while for 0.5 < uj < 1, A converges to a finite value. 
This result suggests that the transition is discontinuous 
(continuous) for < ui < 0.5 (0.5 < uj < 1). 

The case uj = 1: This case is exactly solvable as the 
case of the Smoluchowski coagulation equation Q. We 
consider an arbitrary initial condition of n s (0). In this 
case, c(t) = X^s sn sW ^ s conserved as c(t) = 1, but the 
first moment M\{t) = sn s (t), with the sum exclud- 
ing the largest cluster, is not. The generating function 
gi(/i,i) = J2 S sn s(t) exp(/is) satisfies the relation 



9l = 2( 9l - l)g[, 



(5) 



where the dot (prime) is the derivative with respect to 
time t (/i). Then, g\ is the solution of 



gi(jt,t) = 1 - H(-fi - 2t(gi(ji,t) - 1)), 
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FIG. 4: (Color online) Test of the discontinuous or continuous 
PT for the cluster aggregation network models with various 
uj values. When < uj < 0.5, A decays with increasing iV; 
however, when 0.5 < ui < 1, it converges to a finite value. 
The straight lines are guidelines for eyes. 
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FIG. 5: (Color online) Plot of the giant cluster size G versus 
time for the ER network with different initial conditions. In 
(a), the cluster-size distribution n s (0) is flat with different 
cutoff values s m = N v , with rj = 0, 0.1, 0.2, and 0.3 from right 
to left. In (b), n s (0) decays according to a power law with 



exponent r = 3, 2.5, 
both (a) and (b). 



2, and 1.5 from right to left. JV = 10' in 



where H(fi) = 1 — g\{— /i, 0) is fixed by the initial con- 
ditions of n s (0). The giant cluster size G(t), defined as 
G = 1 — <?i(0 - ,t), can be solved by the self-consistent 
equation G = H(2tG). The obtained G(i) has the form 
near t r as 



G(t) 



2M|(0) 
M 3 (0) 



(2M 2 (0)f - 1) . 



(7) 



where M„(0) = X)' s s ' ln s(0) is the initial n-th mo- 
ment. Also, one finds that the second moment M 2 (t) = 
Y^' s s 2 n s (t), obtained from g[(0~ ,t), behaves as 



M 2 (t) = 



M 2 (0) 



l-2M 2 (0)i| 



(8) 



for t <t c = 1/(2M 2 (0)), and for t > t c as t -> t+. 

These solutions for arbitrary initial conditions are used 
to study the first-order transition in the ER network be- 
low. For the cluster aggregation network model, the ini- 
tial condition is n s (0) = 5 Sll . Then, M 2 (0) = M 3 (0) = 1, 
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TABLE I: When the number of cluster sizes at initial time obeys a power law n a (0) = As~ T for s = 1, . . . ,s„ 
amplitude A, the second moment at initial time M2(0), the critical point t c , and the critical behavior of the giant cluster size 
G(t). Type of PT is specified for each case. The listed t c and G are the ones in the thermodynamic limit. 

T A A/ 2 (0) tc G(t > t c ) type of PT 

(i) < r < 2 4^r (2 7 T)s "' 1 discontinuous 

fii) r = 2 t-^ — , " m 1 discontinuous 

v / In s m In s m 

1 s 3-7 " t-2 

(iii) 2 < r < 3 t? — 777= — r oc t 3 ~ T continuous 

v ' C(t-1) CO-1)(3-t) 

(iv) r = 3 c(t 1 _ 1) oc Jre~ c(2)/2t continuous 

(v) 3<r<4 -^-^ l^iy 2^-2) oc(t-t c )~ continuous 

( vi ) r = 4 TTT^TT 77§T l®T « TnTT^iTT continuous 



2C(2) ln(t-t c ) 
C(r-l) 

C(r-l) C(r-l) 2C(r-2) 



(vii) r > 4 ffr-n 2 C Ar-2) cct-tc continuous 



and consequently, t c = 1/2, which is the well-known ER 
value. The giant cluster size exhibits a continuous tran- 
sition at t c with n B (t c ) ~ s -5 / 2 . 

It is often the case that starting from n s (0) = 5 s .i, 
the cluster-size distribution n s (t) exhibits a power-law 
behavior (or with hump) in s just before or at the tran- 
sition point, even when the dynamics is different from 
the ER (see Fig. [2|. To see how such n s evolves un- 
der the ER dynamics from then on, we consider here 
two particular cases in which M2(0) and M3(0) depend 
on N. First, we assume that n s (0) follows a flat distri- 
bution, n s (0) = no, in the range < s < s m , where 
s m , the size of the largest cluster at t = 0, depends 
on JVass m = N n . Then, n = 2N^, M 2 (0) cx 
JV", and M 3 (0) cx 7V 2r '. Then, a PT takes place at 
t c (N) = 1/2M 2 (0) cx AT-'?, and G(t) - r(2M 2 (0)f - 1) 
for t > t c (N) from Eq. ([7]), where r turns out to be in 
0(1). Thus, if time t is scaled as t' = tM 2 (0), then 
one can show that G(t') is the solution of G = H(3t'G), 
where H(x) = 2 ^ =1 (-l) Tl+1 x ,l /(n!(n + 2)) is a regu- 
lar function qualitatively similar to H(x) = 1 — e~ x of 
the standard ER problem. Hence, G(t') has a mean field 
behavior similar to the original ER case. This scaling be- 
havior implies that while SG(t) = G(ti) — G(to) increases 
by 0(1), A = h - 1 does so by - 0{N^). Thus, we 
have a first-order transition as N — > 00 (Fig.(5Ja)). 

Second, we suppose that the initial condition is given 
as n s (0) = As~ r in the range < s < s m , where A is 
the normalization constant determined by the condition 
y~] sn s = 1. A is given in Table I for various ranges of 
t, together with the initial second moment for arbitrary 



s m , the critical point, and the giant cluster size G(t), 
obtained from Eq. © 

In short, the PT occurs at t = for r < 3, and at finite 
t c for t > 3. This behavior is related to the divergence 
of the second moment M 2 (0) since time is scaled in the 
form t' = 2tM 2 (0). The transition is discontinuous when 
t < 2, but continuous when 2 < t < 4. This difference 
originates from the fact that the ratio M 2 (0)/M3(0) is 
finite for the former, while it vanishes for the latter. For 
t > 4, both M 2 (0) and Ma(0) are finite, resulting in the 
classical percolation behavior at a finite t c . 

In summary, we have introduced a cluster aggregation 
network model, in which discontinuous percolation 
transitions occur when the connection kernel is sub- 
linear as Kij ~ [ij) u with < lj < 1/2. Even for the 
ER network, a discontinuous PT can also be obtained 
by using initial conditions where M 2 (0) diverges and 
M% (0)/Ms (0) remains finite. The simple model mani- 
fests explicitly the role of the abundance of large-size 
clusters just before a transition point as a mechanism of 
the discontinuous PT [l7|. We expect that the cluster 
aggregation network model can be used to study under- 
lying dynamics of the explosive percolation transition 
of random ER network under the Achlioptas process [l3| . 

This work was supported by a KOSEF grant Acceler- 
ation Research (CNRC) (Grant No.Rl 7-2007-073-01001- 
0) and by the NAP of KRCF. 

Note added.-After the submission of this paper, we 
became aware of a work [20j | , which starts from the same 
motivation as ours. 
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